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ABSTRACT 

Context. There is still considerable debate over the progenitors of type la supernovae (SNe la). Likewise, it is not agreed how single 
white dwarfs with masses < 0.5 M Q can be formed in the field, even though they are known to exist. 

Aims. We consider whether single low-mass white dwarfs (LMWDs) could have been formed in binary systems where their com- 
panions have exploded as a SN la. In this model, the observed single LMWDs are the remnants of giant-branch donor stars whose 
envelopes have been stripped off by the supernova explosion. 

Methods. We investigate the likely remnants of SNe la, including the effects of the explosion on the envelope of the donor star. We 
also use evolutionary arguments to examine alternative formation channels for single LMWDs. In addition, we calculate the expected 
kinematics of the potential remnants of SNe la. 

Results. SN la in systems with giant-branch donor stars can naturally explain the production of single LMWDs. It seems difficult 
for any other formation mechanism to account for the observations, especially for those single LMWDs with masses < 0.4 M Q . 
Independent of those results, we find that the kinematics of one potentially useful population containing single LMWDs is consistent 
with our model. Studying remnant white-dwarf kinematics seems to be a promising way to investigate SN la progenitors. 
Conclusions. The existence of single LMWDs appears to constitute evidence for the production of SNe la in binary systems with a 
red-giant donor star. Other single white dwarfs with higher space velocities support a second, probably dominant, population of SN 
la progenitors which contained main-sequence or subgiant donor stars at the time of explosion. The runaway stars LP400-22 and US 
708 suggest the possibility of a third formation channnel for some SNe la in systems where the donor stars are hot subdwarfs. 
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1. Introduction 

Type la supernovae (SNe la) are of major astrophysical impor- 
tance. They have acquired particular cosmological significance 
since they have been used to measure the expansion history of 
the Universe (Riess et al. 1998; Perlmutter et al. 1999; Riess 
et al. 2004). Understanding their nature is also of importance 
for understanding the metallicity evolution and star-formation 
history of galaxies (e.g. Canal, Ruiz-Lapuente & Burkert 1996; 
Matteucci & Recchi 2001). Despite their importance, there is 
still no agreement on the nature of their progenitors. 

There is broad agreement that the destruction of a white 
dwarf (WD) in a thermonuclear explosion constitutes the su- 
pernova event itself, but there are two main classes of com- 
peting models for the events which lead to the explosion. In 
the single-degenerate scenario, the doomed WD accretes mat- 
ter from a non-degenerate companion (Whelan & Iben 1973; 
Nomoto 1982; Han & Podsiadlowski 2004). In the double- 
degenerate scenario, the mass donor is a second WD; the most 
commonly considered scenario involves the merger of two CO 
WDs (Iben & Tutukov 1984; Webbink 1984; also see Martin, 
Tout & Lesaffre 2006 for a variant of this scenario). An ex- 
plosion following the merger of two WDs would leave no rem- 
nant, whilst the companion star in the single-degenerate scenario 
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would survive and be potentially identifiable (Ruiz-Lapuente 
1997; Podsiadlowski 2003; Ruiz-Lapuente et al. 2004). 

There has been no conclusive proof to date that any indi- 
vidual object is the surviving non-degenerate donor from a SN 
la explosion. If Chandrasekhar mass WD-WD mergers do not 
lead to SNe la, they are expected to leave a remnant neutron star 
via acretion -induced collapse (AIC, see Nomoto & Iben 1985; 
Nomoto & Kondo 1991; but see also Yoon, Podsiadlowski & 
Rosswog 2007). At present we do not know whether all WD- 
WD mergers do leave remnants - in which case the double de- 
generate scenario could not be responsible for SNe la - and it 
seems unlikely that that this will become clear in the near future 
(but see, e.g., Levan et al. 2006). Hansen (2003) first noticed that 
observed high-velocity WDs (Oppenheimer et al. 2001) could 
have been produced through SNe la; such WDs would be the de- 
scendants of non-degenerate mass donors in the pre-supernova 
binaries. Hansen's idea seems to be consistent with more de- 
tailed work on the ages of the WDs in the Oppenheimer et al. 
sample by Bergeron et al. (2005) and deserves further atten- 
tion, but by itself it is not a clinching argument for the single- 
degenerate channel. Nor has the evidence that the SN la rate is 
different for different stellar populations (Mannucci et al. 2005) 
led to firm conclusions. The strongest direct evidence that non- 
degenerate donor stars can lead to normal type la supernovae has 
been provided by Patat et al. (2007), who observed circumstel- 
lar material around SN 2006X which seems extremely hard to 
reconcile with a double-degenerate progenitor. 
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Here we suggest that the observed, apparently single, low- 
mass white dwarfs (LMWDs) provide evidence that at least 
some SN la explosions have occurred with non-degenerate 
donor stars. We define LMWDs as WDs which are too low 
in mass to have been produced by single-star evolution as we 
currently understand it. A population of single LMWDs has 
been implied by, e.g., the work of Maxted, Marsh & Moran 
(2000)Q We also investigate the apparently single ultra-cool 
white dwarfs (UCWDs) as potentially containing a useful subset 
of the LMWD population and indicate how further observations 
of the kinematics of this and other populations could lead to con- 
straints on the progenitors of SNe la. 

In sections [2] and [3] we argue that the existence of single 
LMWDs is most naturally explained by the single-degenerate 
model for SNe la. In section |4] we introduce UCWDs, discuss to 
what extent the observed single UCWDs might be a useful sam- 
ple of LMWDs, and consider in what way the observed single 
UCWD population is consistent with a SN la origin. 

2. Formation channels for single LMWDs 

Current evidence suggests that 1 M zero-age main-sequence 
(ZAMS) stars, left to evolve in isolation, produce white dwarfs ^ 
0.55 M (e.g. Han, Podsiadlowski & Eggleton 1994; Weidemann 
2000). In order to produce low-mass helium WDs (with masses 
< 0.5 M ), it is necessary to remove the envelope of a star be- 
fore it is able to ignite helium. Within the age of the Universe, it 
is almost certainly impossible for a single star to produce WDs 
with masses close to 0.3 M , or even 0.17 M , as recently in- 
ferred for the runaway WD LP 400-22 (Kawka et al. 2006), or 
0.23 M , the estimated mass of the apparently single UCWD 
LHS 3250 (Bergeron & Legget 2002). 

2.1. Production of single LMWDs by single stars? 

It would be simplistic to conclude from the existence of cur- 
rently single LMWDs that single stars can produce LMWDs. 
However, arguments have been made in favour of a single star 
channel for the production of some LMWDs, at least in signifi- 
cantly metal-rich populations (Kalirai et al. 2007; Kilic, Stanek 
& Pinsonneault 2007; see also Han et al. 1994). Whilst this is 
a possibility, that interpretation is not the only one which can 
explain the observations: this paper provides one alternative ap- 
plicable to both field stars and clusters, and dynamical interac- 
tions in dense clusters may provide another route (e.g. Adams, 
Davies & Sills 2004)0 Bedin et al. (2008) and Van Loon, Boyer 
& McDonald (2008) do not support the single-star formation 
channel for LMWDs suggested by Kalirai et al. (2007). 

Even if some WDs less massive than ~ 0.55 M can be pro- 
duced by single-star evolution with super-solar metallicity, the 
masses are unlikely to approach those of LP 400-22 or LHS 
3250 (as above). Follow-up calculations to Han et al. (1994) im- 
ply that only LMWDs with M > 0.4 M might be produced 
from such a single-star channel, even at high metallicity (Meng, 



1 See also Marsh, Dhillon & Duck (1995); their study was originally 
motivated by their interest in finding double degenerate binaries as po- 
tentional SN la progenitors. 

2 Note that mass segregation tends to move low-mass objects out- 
wards in such clusters, so LMWDs might be formed in the core and 
observed far from the centre. 

3 Bedin et al. argue that there may not, after all, be an unusual single 
LMWD population in NGC 6791 and Van Loon et al. see no evidence 
for enhanced stellar mass loss in infra-red observations of the cluster. 



Chen & Han, 2007). Kilic et al. (2007b) also state that their pro- 
posed single star channel is highly unlikely to be relevant for 
WD masses as low as 0.2 M . 

Nothing of what follows would be significantly affected if 
the maximum mass of the LMWD category had to be revised 
downwards slightly to account for non-standard single-star evo- 
lution. We will mostly be considering LMWDs with masses 
< 0.3 M , at which point any presently proposed single-star 
channel is not expected to contribute. 

2.2. Production of single LMWDs by binary stars? 

In binary systems we can invoke mass transfer (and sometimes 
ablation by a pulsar companion) in order to explain the observed 
binary LMWDs (e.g. van Kerkwijk et al. 2000; Liebert at al. 
2004). Apparently single LMWDs must also be formed within 
an interacting binary system - either we have not detected their 
companion or the binary has been disrupted. One attractive for- 
mation channel stands out: the formation of the LMWD in a bi- 
nary where the binary companion exploded in a SN la. Before 
examining this channel in more detail, we discuss possible alter- 
native explanations for single LMWDs. 

2.2.1. Alternatives to SN la: Core-collapse supernovae? 

The natural alternative to a SN la in explaining the disruption 
of a binary is a core-collapse supernova (chosen such that the 
system becomes unbound). However, forming a LMWD in such 
a system is challenging. The binary must remain intact long 
enough for the WD progenitor to lose its envelope such that it 
will later form a LMWD; this implies that the initially more 
massive star has transferred its envelope to the secondary and 
then becomes a WD. In order to form a LMWD, the primary 
mass must be < 4.5 M , which our own stellar calculations 
find to have a core mass at the end of the main sequence of 
« 0.52 M Q In order to produce a core-collapse supernova (re- 
quiring Mzams S 8 M ), the initially less massive star would 
thus have to accrete the large majority of the envelope of the pri- 
mary soon after the primary has left the main sequence, and the 
initial mass ratio would have to be close to 1 . This highly opti- 
mistic scenario does not produce a distinctly low-mass WD; to 
produce a 0.3 M WD in this way requires a ZAMS mass for the 
primary of w 3 M , precluding a core-collapse supernova in the 
binary system. 

Our arguments above are generous; Davies, Ritter & King 
(2002) found a much more restrictive result. In a different con- 
text they investigated the evolution of systems where a WD is 
formed before a core-collapse supernova occurs in the system. 
They concluded that the mass of the white dwarfs generated 
in this way, Mwd k. 1 M '. 

2.2.2. Alternatives to SN la: Circumbinary discs? 

If cataclysmic-variable (CV) evolution is driven by circumbinary 
discs, the donor star may eventually be entirely consumed (the 
'White Widow' scenario; see Spruit & Taam 2003, following 
Spruit & Taam 2001 and Taam & Spruit 2001). For this mech- 
anism to explain single LMWDs, the WD in the progenitor CV 
must also have been a LMWD; the WD may not gain much mass, 



4 We use Eggleton's stellar evolution code (Eggleton 1971; Pols et al. 
1995) with a metallicity of 0.02 along with the convective overshooting 
calibration of Pols et al. (1998). 
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as the matter it accretes can be ejected via nova explosions, but 
it is unlikely to become significantly less massive. 

It is not clear whether this mechanism operates in CVs: at- 
tempts to detect circumbinary discs have inferred disc masses 
several orders of magnitude below the required values (Dubus 
et al. 2004; Muno & Mauerhan 2006). Hence in the abscence of 
further supporting evidence we consider this potential formation 
channel unlikely at present. 

2.2.3. Alternatives to SN la: Accretion-induced collapse? 

Some systems will contain an accreting WD which succeeds in 
reaching the Chandrasekhar mass but fails to produce a super- 
nova as the WD is predominantly composed of oxygen, neon & 
magnesium (ONeMg) rather than carbon and oxygen (CO). This 
can occur either because the WD began accreting as an ONeMg 
WD or because the accretion rate onto the WD did not allow the 
WD to remain a CO WD (e.g. Nomoto & Iben 1985; Nomoto 
& Kondo 1991; Martin et al. 2006). Such WDs will produce a 
neutron star (NS) via AIC. Currently it does not appear likely 
that AIC produces sufficiently large kicks to disrupt such close 
binaries (see, e.g., Podsiadlowski et al. 2004). 

2.2.4. Alternatives to SN la: White-dwarf mergers? 

Single LMWDs may be the product of the merger of two low- 
mass He WDs, with formation rates comparable to or greater 
than the SN la rate (see, e.g., Han et al. [2002] and references 
therein). However, Han et al. predict masses in excess of 0.4 M . 

3. Single-Degenerate SN la Populations 

3.1. Expected formation channels 

We do not present an exhaustive description of the full evolu- 
tionary histories for single-denegerate SN la progenitors (see, 
e.g. Whelan & Iben 1973; Nomoto 1982; van den Heuvel et al. 
1992; Rappaport, Di Stefano & Smith 1994; Hachisu, Kato & 
Nomoto 1996, 1999; Li & van den Heuvel 1997; Langer et al. 
2000; Hachisu & Kato 2001 ; Han & Podsiadlowski 2004). There 
is no clear consensus on which donor stars are likely to produce 
a type la supernova. The favoured options involve either donors 
on the main sequence (MS) or the subgiant branch (known as the 
supersoft channel), or red-giant (RG) donors. 

While the supersoft channel (e.g. Han & Podsiadlowski 
2004) is arguably the favoured channel for the majority of SNe 
la, Hachisu, Kato & Nomoto (1996, 1999) and Hachisu & Kato 
(2001) suggest situations in which a low-mass giant star may 
take a WD to the Chandrasekhar mass Mch at long orbital pe- 
riodsQ Sokoloski et al. (2006) used the 2006 outburst of RS 
Ophiuchi to confirm the conclusions of Hachisu & Kato by infer- 
ring that RS Oph contains a very massive WD (Mwd ~ 1 -4 M ). 
It is worth noting that we cannot be sure that RS Oph contains 
a CO WD rather than an ONeMg one and so we cannot be sure 
that it will explode rather than collapse. Observational support 
for a giant donor in a system which produced a SN la has been 
provided via the observations by Patat et al. (2007) of SN 2006X. 

King, Rolfe & Schenker (2003) have also suggested that 
an accreting WD may not reach Met, via the supersoft chan- 
nel alone, but that a later phase of WD growth could occur in 

5 The point at which an explosive nuclear runaway occurs in a non- 
rotating CO WD is slightly below the Chandrasekhar mass: Nomoto, 
Thielemann & Yokoi (1984) calculated a mass of ~ 1.378M G . 



long-period dwarf novae. They argue that, even though the av- 
erage mass-transfer rate does not reach the steady-burning band 
(Paczyriski & Zytkow 1978; Nomoto & Kondo 1991), the ac- 
cretion rate may be high enough for the WD to grow during 
dwarf nova outbursts driven by the thermal-viscous disc insta- 
bility (Cannizzo, Ghosh & Wheeler 1982). Providing the correct 
mass-accretion rate for the CO WD to grow to Men is a signifi- 
cant uncertainty in all these models. 

3.2. Remnant mass 

In order to understand the formation of LMWDs in systems 
which produce SN la explosions, we must consider the mass and 
evolutionary stage of the donor star at the point of the explosion 
and also the extent to which the donor loses mass because of the 
explosion. There is a clear division between pre-giant and giant 
donor stars, with giant donors apparently able to leave LMWD 
remnants. 

Marietta, Burrows & Fryxell (2000) performed numerical 
simulations of the effect of a SN la explosion on the compan- 
ion star. They found that 0.15 to 0.17 M Q is stripped away 
from a 1 M main-sequence or subgiant companion by the high- 
velocity ejecta. Han & Podsiadlowski (2004) found in their pop- 
ulation synthesis simulations of the supersoft channel that, at 
the time of the explosion, the companion has a mass between 
> 0.5 M and 2.2 M , with a typical mass of 1 M (for more 
details see also Han 2008). Applying the results of Marietta et 
al. as a percentage - 15% of the donor mass - leads to a low- 
est estimated remnant mass of « 0.42 M . If the WD explodes 
as it reaches Mch, then this remant mass is a lower limit for 
the MS channel, assuming negligible subsequent mass loss in 
a wind. Hachisu & Kato (2001) found a lower limit on the mass 
of the donor from the supersoft channel (at the time of the SN) 
of > 1.3 M (assuming an initial white dwarf mass of 1 M ). 
Despite these differences, both studies suggest that it is difficult 
to produce LMWDs via main-sequence or subgiant donors. 

Marietta et al. also found that a red-giant donor will lose 
almost its entire envelope (96% - 98%) due to the impact of 
the SN la explosion and leave only the core of the star, provid- 
ing a possible pathway for the formation of a subset of single, 
low-mass He WDs. For the RG channel, Hachisu & Kato (2001) 
found a lower limit on the total donor mass of > 0.4 M . If the 
RG channel produces SNe la, then ram-pressure stripping of the 
donor's envelope would be expected to lead to the formation of 
LMWDs. The remnant WD mass is dependent on the core mass 
of the donor at explosion and is therefore strongly correlated 
with the orbital period (see section [3~4l . 

One formation channel that is rarely discussed in the litera- 
ture is one where the donor star is a hot subdwarf star (see, e.g., 
Geier et al. 2007). We do not expect significant stripping of the 
donor by the supernova ejecta in this case, as the donor star will 
be tightly bound, but the mass of the donor star could easily be 
low enough for a single LMWD to be formed by the natural evo- 
lution of the donor star. 

3.3. Sub-Chandrasekhar mass explosions 

A variation on the above models for SN la involves the explo- 
sion of sub-Mch CO WDs covered with a thick helium layer 
(Woosley & Weaver 1994). In that model, the detonation of the 
helium layer is responsible for triggering the supernova. Fink, 
Hillebrandt & Ropke (2007) found that sub-Mch explosions 
were unlikely to be able to explain either normal or subluminous 
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SN la, but there could be implications for our LMWD formation 
channel if a significant fraction of SNe la were found to be pro- 
duced by sub-Mch detonations^ Qualitatively there would be lit- 
tle change to our model, as the RG donors would still be stripped 
of their envelopes and produce LMWDs. The quantitative rem- 
nant mass distribution may be different. For example, these sub- 
Mch explosions might plausibly happen when the donor stars 
are lower on the giant branch than for the standard model. In 
that case, the typical remnant WDs may be less massive, the or- 
bital periods at explosion lower and the runaway remnant veloc- 
ities higher than for detonations at Men- The remnant velocity 
is considered hereafter as a diagnostic of the orbital period at 
explosion, assuming Chandrasekhar-mass explosions. 

3.4. Binary orbits and runaway velocities 

If single LMWDs have been released from binary systems in 
which the other component has exploded as a type la supernova, 
the space velocity of the remnant should be a useful diagnostic of 
the orbital period at explosion. As our arguments above suggest 
that the LMWDs are most likely to originate in systems with red- 
giant donor stars, the relationship between core mass and orbital 
period in such systems can act as a further constraint. 

In what follows we assume that the donor stars are filling 
their Roche lobes, as it seems to us that the mass transfer in sys- 
tems which produce SN la is most likely to be due to Roche-lobe 
overflow. However, it is not known whether this is the case as, 
for example, it is unclear whether the donor star in RS Ophiuchi 
is rotating synchronously with the binary orbit (see, e.g. Murset 
& Schmid 1999; Zamanov et al. 2007). If the donor stars do not 
fill their Roche lobes, then the following method may slightly 
overestimate the runaway velocities for a given remnant mass. 

If we define q as MijM\ (where M\ is the mass of the SN 
la progenitor and Mi is the mass of the companion producing 
the LMWD), write the total mass of the system (in solar units) 
as M tot , and the pre-SN orbital period (in days) as Pdays, we can 
write the pre-SN orbital velocity as: 



213 / M, 



orb 



1 + q \ Pi, 



1/3 



km/s. 



(1) 



This shows that, for M2 <K Mi , the orbital velocity is relatively 
insensitive to the donor mass. 

Furthermore, assuming that the remnant WD mass equates to 
the core mass of the donor at the time of the SN explosion, then 
the well-defined relationship between core mass and radius for 
giant-branch donors leads to an expression for the orbital period 
of the system at the supernova stage. Rappaport et al. (1995) 
found the period-mass relation: 



0.374 



R Q M, 



.5 

wd 



1 +4M\ 

wd 



+ 0.5 



3/2 



M~ l J 2 days, 



(2) 



where M w( j is the mass of the future WD (currently the core of 
the giant star) in units of solar masses. Their preferred value for 
the fitting parameter Ro was 4950 R Q , which we also adopt. If 
we now assume a primary mass of 1 .4 M and a total secondary 
mass at the time of explosion of 0.5 M Q , we obtain for a range 
of core masses: 



6 If such sub-Mch explosions occur but are not seen as SN la, it is 
not clear to us whether such an event would remove the RG donor's 
envelope in order to produce a LMWD. 



P orb (0.2 Mo) = 6.74 d 

P orb (O.25M ) = 23.7d 
P orb (0.3 Me) = 69.4 d 
P OTb (O.35M ) = 171.8 d 



Vorb = 103 km s 
V orb = 68 km s" 1 
y orb = 47 km s" 1 
V orb = 35 km s" 1 



where V orb is the total orbital velocity. Using the same com- 
ponent masses, we record for comparison that P orb = 1 d cor- 
responds to V OI b ~ 194kms~ 1 and f orb = 1 h to V olb ss 
560 km s . The orbital velocities at explosion in the simulations 
of Han & Podsiadlowski (2004) range from 80-230 km s for 
MS and subgiant donors. 

RS Ophiuchi has an orbital period of »457 d. Inverting equa- 
tion [2] above, this corresponds to a core mass for the donor of 
slightly over 0.4 M , still within the mass range for a LMWD 
should the envelope be removed. 

The arguments in section [372] suggest that LMWDs are pro- 
duced by giant donors, as long as the orbital period is not so 
long that the core has already grown to 0.5 M Q by the time their 
tenuous envelopes are stripped by the supernovae ejecta. So, in 
contrast to the high-velocity WDs observed by Oppenheimer et 
al. (2001) and interpreted by Hansen (2003) as remnants of SNe 
la with main-sequence donors, it would be consistent to find 
that the single LMWD population was not significantly kinemat- 
ically heated. 

In section 13.21 we identified two potential SN la formation 
channels able to produce LMWDs. We note that whilst one of 
those sets of donor stars (red giants) would leave LMWD rem- 
nants with low runaway velocities, the other (hot subdwarfs) 
would result in high-velocity LMWDs. 

4. Single UCWDs as LMWDs & SN la remnants? 

We have argued that single LMWDs can be produced from 
single-degenerate SNe la with red-giant donors. Single LMWDs 
are inferred to exist and the most natural explanation, espe- 
cially for the lower-mass LMWDs, seems to be that some single- 
degenerate SN la occur with red giant donors. However, there is 
no obvious collected sample of LMWDs to examine as potential 
SN la remnants. Independently of our arguments above, it may 
well be that the known set of apparently single UCWDs con- 
stitutes or contains a useful sample of single LMWDs. When a 
suitable sample of single LMWDs becomes available our work 
should be extended. 



4.1. The UCWD sample 

In selecting sub-samples of the WD population, UCWDs (see, 
e.g. Harris et al. 1999 & 2001; Gates et al. 2004; Wolf 2005) are 
clear outliers in a colour-colour diagram (see Fig . [TJ . Their opti- 
cal colours distinctly separate them from the normal WD popula- 
tion, and they are an interesting curiosity in appearing to become 
bluer as they cool down, possibly due to the effects of collision- 
ally induced absorption (CIA) by hydrogen molecules in the at- 
mosphere (Bergeron et al. 1994; see also Kowalski & Saumon 
2006). When CIA affects only the infrared part of the spectrum, 
WDs are classed as cool; if CIA also affects the optical colours, 
then the WD is admitted into the select group of UCWDs (see 
Fig. 1 ; also Wolf 2005). The transition temperature between cool 
and ultra-cool is <; 4000K. 

The known UCWDs constitute a clean observational sample; 
they are easy to identify, and with such low luminosities (no- 
tably in the ultra-violet and near infra-red) and line-free spectra, 
it would be hard to hide a light-emitting close companion that is 
anything but another UCWD (see section l4~3l . Table Q] contains 



Stephen Justham et al.: On type la supernovae and the formation of single low-mass white dwarfs 



5 



Table 1. Distances and tangential velocities of the UCWD sample (as in Fig. [TJ). For all objects except LHS 3250, the distance 
is inferred from an assumed absolute magnitude. Our preferred measure assumes LHS 3250 is a typical UCWD; we also present 
alternative values based on the range of absolute magnitudes assumed in the literature. 



Name 



LHS 3250 1 
LHS 1402 2 
SDSS J0947 3 
SDSS J1001 3 
SDSS J1220 3 
SDSS J1337 4 
SDSS J 1403 3 
COMBO-17 J1143 



Estimated properties for My = 15.7 " 
Distance (pc) Tangential velocity (km/s) 



Estimated properties for My = 16. 5 +1.0 " 
Distance (pc) Tangential velocity (km/s) 



30 

24 c 

47 

64 

64 

54 

44 

169 



81 
56 
18 
107 
154 
46 
60 
42 



LHS 3250 has a known parallax 
12-31 d 28-72 
21-52 8-21 

28-71 47-119 

28-71 68-170 
23-58 20-49 
19^19 26-66 
37-188 9^6 



": i.e. taking LHS 3250, with a known trigonometric parallax, as representative of all UCWDs. 

*: i.e. using a conservative assumption adopted by Salim et al. (2004), Gates et al. (2004) and Wolf (2005). 

c : For LHS 1402, we compare a B magnitude of 18.32 (Oppenheimer et al. 2001) with LHS 3250 (18.85; Harris et al. 1999). 

d : Correcting for different photometry, we compare a B magnitude of 18.2 with the g magnitude (20.04) of SDSS 1001. 

1 : Harris et al. 1999 2 : Oppenheimer et al. 2001 3 : Gates et al. 2004 4 : Harris et al. 2001 5 : Wolf 2005 




Fig. 1. A colour-colour diagram of white dwarfs. The known 
ultra-cool white dwarfs are represented by large circles, while 
the Kleinman et al. (2004) sample of normal WDs are shown as 
small circles. Table [T]lists the UCWDs plotted here along with 
the relevant references. (Adapted from Wolf [2005].) 



the UCWD sample we use. The estimates of the tangential ve- 
locities for these objects depend on their assumed distances, and 
Table Q] shows the range of velocities obtained for assumptions 
taken from the literature. For that sample we shall adopt the tan- 
gential velocities obtained by taking the absolute magnitude of 
LHS 3250 (which is the only UCWD with a known parallax) to 
be representative of the whole sample. 

Since we first began this work, a new sample of twenty-four 
UCWD candidates has been presented by Vidrih et al. (2007). 
We are not convinced that they are cold enough to conform to 
our strict criteria as UCWDs, which may mean that this sam- 
ple is more likely to be contaminated by non-LMWD objects 
(see section 4.2). However, in our later figures we shall show for 
comparison this new, independent, sample alongside the smaller 
set from Table [TQ 



4.2. Are the observed UCWDs mostly LMWDs? 

The one UCWD with a known parallax (LHS 3250) has an abso- 
lute magnitude of My = 15.72. This is brighterthan expectedfor 
anything other than a LMWD; hence Bergeron & Leggett (2002) 
conclude that the mass of LHS 3250 is 0.23 M o Unfortunately 
we do not have such good mass estimates for all single UCWDs. 

The argument that the observed UCWDs are mostly 
LMWDs is partly built upon theoretical WD cooling tracks. 
Whereas a 0.6 M Q WD takes more than 9 Gyr to cool to 4000 K 
{after the formation of the WD), the 0.3 M WD of equivalent 
composition takes less than 4 Gyr (Bergeron, Leggett & Ruiz 
2001; Bergeron et al. 2005). As the cooling of WDs is a function 
of composition it is likely that not all UCWDs are LMWDs (see, 
e.g., Hansen 1999, who requires WD masses < 0.25 M in or- 
der for those objects to cool to 4000K within 7 Gyr) Extremely 
low-mass WDs (< O.17M ) seem to cool more slowly than more 
massive WDs, due to the retention of a relatively thick hydro- 
gen envelope (e.g. Panei et al. 2007). It is not clear whether 
the LMWD remnants that have been formed by having their en- 
velopes forcibly removed by a supernova Shockwave will retain 
a thick hydrogen envelope. 

Hence we expect that the observed, apparently single, 
UCWDs are dominated by single LMWDs if they exist, partly 
as they are significantly more luminous than massive WDs and 
hence more likely to be discovered. Although it is unfortu- 
nate that we cannot prove what fraction of single UCWDs are 
LMWDs, in section l4~4l we show that the observed numbers of 
single UCWDs could all be single LMWDs produced via a SN 
la explosion. 



4.3. Are the apparently single UCWDs mostly single? 



Seven of the eight UCWDs in Fig. 
ion. The exception - SDSS J0947 - 



1 have no known compan- 
has a common proper mo- 



7 We have not included an extremely recent sample of seven very 
cool white dwarfs (Harris et al. 2008). 



8 Harris et al. (2001) also suggest that LHS 3250 is probably a low- 
mass helium core WD based on its absolute magnitude; however, the 
earlier paper by Harris et al. (1999) noted that a pair of UCWDs in a 
binary might help explain the higher luminosity without recourse to a 
low-mass WD. 

9 We note in advance that this mass of < 0.25 M G is in good agree- 
ment with our results in section 14.61 combined with the core-mass 
orbital-period relation; it is thus self-consistent. 
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tion companion (Gates et al. 2004). The 20 arcsecond angular 
separation of J0947 from its potential companion implies a pro- 
jected separation of over 2 x 10 5 R for a distance of 47 pc; if 
that companion really forms a binary with J0947 (rather than be- 
ing chance projection), it could not have influenced the evolution 
of the progenitor of the J0947 UCWD. 

We do, however, need to consider whether these apparently 
single UCWDs really are single. We fully expect that UCWDs 
should exist in binary systems, but in the following we argue that 
these are unlikely to contaminate our sample; in many cases, a 
companion would even completely hide a UCWD. 

4.3.1. Non-degenerate companions 

An M-dwarf with an absolute V magnitude of ~ 16 - similar 
to the UCWD LHS 3250 - would have a mass of ~ 0.1 M 
(Delfosse et al. 2000). Such a star would be bright in the infra- 
red, where the emission of UCWDs is strongly suppressed[3 
Furthermore, M-dwarfs are rich in spectral lines, so could be 
easily detected. 

4.3.2. White-dwarf companions 

A binary containing a non-ultra-cool WD should be identifiable: 
an advantage of UCWDs is that their low luminosity and fea- 
tureless spectrum makes it hard for them to possess an unde- 
tectable hotter WD companion. Their characteristic spectral en- 
ergy distributions mean that a non-UCWD companion would be 
brighter than the UCWD in either near infra-red or ultra-violet 
light. However, a spectroscopic UCWD-UCWD binary would be 
difficult to distinguish from a single UCWD. We have no reason 
for thinking that such binary UCWDs do not exist, and Harris et 
al. (2008) have discovered one system which might eventually 
be expected to become such a binary. Although we do not ex- 
pect that they are common enough to dominate the population, 
the possibility that the apparently single UCWDs have extremely 
cool or faint WD companions should be studied further. 

4.3.3. Neutron-star companions 

From the arguments in section 2.1 it is clear that, in order to 
make a LMWD in a system where the NS was produced in a 
core-collapse supernova, the LMWD progenitor must lose mass 
after the formation of the NS. Hence, either mass transfer onto 
the NS or a common-envelope phase (Paczyhski 1976) would be 
required. In the former case, the pulsar is expected to be recycled 
into a millisecond pulsar (MSP). 

It is thus reasonable to expect that UCWDs might be seen 
as companions to millisecond pulsars. This has been claimed 
for PSR J075 1 + 1807 (Bassa, van Kerkwijk & Kulkarni 2006a), 
although the LMWD in that system is not quite formally ultra- 
cool. The nearby millisecond pulsar J0437-47 1 5 also seems to be 
accompanied by a borderline UCWD (Danziger, Baade & Delia 
Valle 1993). It could be argued that, since we do see LMWDs 
with pulsars, it is not a surprise that some of these LMWDs hap- 
pen to be in systems where we do not see the pulsar (see also van 
Kerkwijk, Bergeron & Kulkarni 1996; van Kerkwijk et al. 2000; 
Bassa et al. 2006b). However, the space densities of these dif- 
ferent samples are extremely different. The companions to pul- 



10 For example, a ~ 0.1 M G star is about six magnitudes brighter in 
the J band than the V band (Delfosse et al. 2000), wheras LHS 3250 
and SDSS 1337 are dimmer in the J than V bands (Harris et al. 2001). 
This trend is greater at longer wavelengths. 



sars have been discovered because they are orbiting a pulsar, 
unlike the photometrically-discovered UCWD sample: Bassa, 
van Kerkwijk & Kulkarni state that PSR J075 1 + 1807 is -0.6 
kpc away, considerably further than any of our UCWD sample 
(see Table [B, and PSR J0437-4715 is 150 pc away (Danziger et 
al. 1993), beyond all but one of the photometrically-discovered 
UCWDs. 

A crude comparison between the possible detection vol- 
umes of the single UCWDs and the pulsar companions gives 
~ (600/70) 3 ~ 630. Here we have considered that COMBO- 
17 Jl 143, as a serendipitously discovered object, does not give a 
good measure of the systematic detection volume of the UCWDs 
in our sample. This suggests that, even if the companions to 
PSR J075 1 + 1807 and PSR J0437-4715 were cool enough to be 
UCWDs, the seven systematically discovered UCWDs in our 
sample are ~ 2200 times more abundant than those with pul- 
sar companions. This factor can be increased by another order 
of magnitude if one assumes that all the known nearby pulsars 
have been studied, but that only a tenth of the sky has been sur- 
veyed to the same depth as the field that produced the SDSS 
UCWD sample. To only observe one pulsar in ~ 20, 000 suppos- 
edly MSP-containing objects would imply extremely narrow- 
beam pulsar emission. Unless those unseen MSPs were some- 
how unusual, this would suggest that many more pulsars ex- 
ist than we currently expect, significantly worsening any mis- 
match between the inferred birthrates of LMXBs and MSPs (e.g. 
Kulkarni & Narayan 1988; Lorimer 1995; Pfahl, Rappaport & 
Podsiadlowski 2003)0 

If a NS was formed in the system through AIC and was not 
subsequently spun up, it could reasonably be expected not to 
emit pulsar radiation. A black-hole companion would, of course, 
not be expected to emit pulsar radiation, but a local space density 
for such black-hole binaries of ~ 10~ 5 pc" 3 (see section |4~4| | is 
highly unexpected (see, e.g., Romani 1998). 

4.4. UCWD population numbers 

Gates et al. (2004) estimated a space density for UCWDs of 3 x 
10~ 5 pc 3 from a sample of 6 objects found in the Sloan digital 
sky survey[3 This rough figure does compare with an estimate 
of the SN la rate integrated over time and space. The local stellar 
density of 0.1 M Q pc~ 3 (Binney & Memfield 1998), combined 
with a mass for the thin disc of ~ 4x 10 1() M and the assumption 
that the mass fraction of UCWDs is constant throughout the disc, 
produces an estimate for the number of UCWDs of ~ 10 7 . 

We approximate the current SN la rate in the Galactic disc 
as ~ 7 x 10~ 3 /yr (using the same disc mass as above and the SN 
la rate per unit mass of Mannucci et al. 2005). Multiplying this 
rate by a Galactic age of ~ 10 Gyr leads to an estimate of a total 
of ~ 7 x 10 7 remnants Q 



11 If we repeat this estimate for the new Vidrih et al. sample of 24 
UCWDs, which extends to a distance of 180 pc but over only 250 square 
degrees of the sky, then we have ~ (600/180) 3 x (24/2) x (41523/250), 
i.e. single UCWDs outnumber those with MSPs by a factor of over 
70,000 (over 27,000 if UCWDs are twice as bright as Vidrih et al. as- 
sume and their survey depth extends to 250 pc). 

12 For the Viridh et al. sample, we estimate a space density of ~ 16 x 
10~ 5 pc~ 3 for a survey depth of 180 pc, and ~ 6 x 10~ 5 pc~ 3 for a survey 
depth of 250 pc. 

13 The supernova rate may well have been different in the past. Hansen 
(2003) made estimates for the number of SN la remnants in the Milky 
Way based on the amount of iron in the Galaxy and found a range be- 
tween 4 x 10 7 and 2.2 x 10 8 objects. 
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Fig. 2. Comparison of the tangential velocities resulting from our Galactic integrations (smooth black curves), with the tangential 
velocities of the observed UCWDs (grey step functions). The solid light grey step function uses the distance estimates from the 
assumption that LHS 3250 is representative of the UCWD population, and the broken grey step functions encompass the wider 
range of distances in Table Q] The dark grey step function uses the new sample of UCWD candidates from Vidrih et al. (2007), 
adopting their assumptions for the UCWD distances. The solid black curves in each panel are for orbital periods at the time of the 
explosion of 100, 10 and 1 days (left to right). The dashed black curve represents a population which receives no 'kick' due to the 
break-up of a binary system. The Kolmogorov-Smirnov acceptance probabilities for the individual models are given in each panel, 
compared to both the objects in Table Q] and the Vidrih et al. sample. Here refers to the model with no added 'kick' velocity 
from binary break-up and Pino, Pio and Pi to the curves representing the 100, 10 and 1 day orbital-period populations, respectively 
(see text). The simulations assumed initial disc scale-heights of 75 pc, 500 pc and 4 kpc (left to right). If these objects were once 
the cores of red-giant donor stars in SN la producing systems, the simplest expectation would be for orbital periods of ~100 d (see 
text). 
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Fig. 3. As Fig. [2] for simulations that do not include any scattering from giant molecular clouds. The plots assume initial disc 
scale-heights of 75 pc, 200 pc and 1000 pc (left to right). 



Even if we halve this number of remnants to allow for the 
WD cooling time (see section l4T2l ). the number of single UCWDs 
is easily consistent with them being produced through SN la. 
Indeed, these numbers suggest that only a subset of SNe la pro- 
duces single UCWDs, which is as expected if only a subset of 
the SN la formation channels can produce LMWDs (see sections 
|3~T1&|121>. 

4.5. Population kinematics: Method 

The lack of lines in UCWD spectra means that we do not know 
the radial velocities for our sample. However, we can exam- 



ine the population kinematics using only the information from 
the tangential velocities. We now investigate whether their ob- 
served space velocities are consistent with single LMWDs re- 
leased from binaries with a range of orbital periods. 



For a range of inital parameters, we integrated the mo- 
tion of 10 5 assumed SN la remnants for up to 10 Gyr through 
the Galactic potential (using a similar procedure to Brandt & 
Podsiadlowski 1995), orientating the orbital velocity vector of 
the donor at random when the binary is disrupted. For each inte- 
grated population, we used a single value of orbital period (and 
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hence orbital velocity) at the time of the explosion^ For calcu- 
lating the orbital velocity at a given orbital period, donor stars are 
assumed to be 0.5 M at the time of the explosion, and the WDs 
are assumed to explode at a mass of 1.4 M . Equation Q] shows 
that our results should be relatively insensitive to those assump- 
tions, but in the future we intend to perform this procedure using 
the output of our binary population synthesis calculations. 

Each remnant is initially located at random within an ax- 
isymmetric Galaxy modelled by two exponential scale-heights 
(vertical and radial). The axisymmetry is also exploited for com- 
putational efficiency: at each integration time-step the view from 
Earth is calculated at all points on the solar circle. A further as- 
sumption is that the remnants can be observed to a distance of 
160 pc - broadly appropriate for UCWDs. Within such a small 
volume, the space velocities should only be a very weak function 
of distance. 

The Galactic potential was taken from Paczyhski (1990), us- 
ing the parameters in Brandt & Podsiadlowski (1995), with the 
addition of scattering from giant molecular clouds (GMCs) ran- 
domly distibuted within the Galactic disc. The total mass in these 
GMCs is assumed to be 1.2 x 10 9 M , with a mass spectrum ex- 
ponent such that diV/dlogm m (see, e.g., Digel et al. 1996; 
Binney & Merrifield 1998). 

We do not include any kick imparted by the supernova ejecta 
in our simulations, since the simulation of Marietta et al. (2000) 
show that the kicks due to the supernova interaction (86 km s _I 
and 49 km s _1 for their main-sequence and subgiant donors, re- 
spectively, with no kick given to the core of the giant donor) 
are generally small compared the orbital-velocity kick from the 
break-up of the binary. 

Since the initial disc scale-height of the progenitor popula- 
tion is uncertain, we present our results for a range of values. 
The vertical scale-height of massive stars is 75 pc (van der Kruit 
1987), and ~ 200 pc is an approximation to a more generic thin 
disc population (e.g., Ojha et al. 1996, who find a scaleheight of 
260 ± 50 pc for the Galactic thin disc and 760 + 50 for the thick 
disc; see also Kroupa, Tout & Gilmore 1993). Given previous 
speculations that UCWDs are so cool because they are very old 
objects, we also modelled initial scale-heights of 500 pc, 1 kpc 
and 4 kpc{3 

4.6. Population kinematics: Results 

Figures [2] &E] present the results of our integrations for orbital 
periods at the time of the supernova of 1, 10 and 100 days, as 
well as a population which received no kick. The population 
which received no kick seems to be difficult to reconcile with the 
observed tangential velocities, except for the most extreme range 
of luminosities consistent with the literature combined with an 
inital vertical scale height of 4 kpc (or greater). The population 
released from a one day orbital period also appears inconsistent 
with the data, whichever initial scale-height is assumed. These 
conclusions are supported by applying the Kolmogorov-Smimov 
test to compare our simulations with the tangential velocity dis- 
tribution produced by assuming LHS 3250 is a typical UCWD 



14 We expect that in reality SN la explosions will occur over a range 
of orbital periods, but the current data only allow us to investigate broad 
trends. 

15 Note that the collective tangential velocities of the sample in Table 
1 are clearly inappropriate for objects from a halo population. 



(TableQ}. Each plot contains a Kolmogorov-Smirnov acceptance 
probability (P x ) for each curve 

Given the considerable uncertainties, we consider that the 
extremely favourable Kolmogorov-Smirnov test for the popu- 
lation released from 100 day orbits with a 4 kpc scale height 
should not be over-interpreted, especially as a real SN la popu- 
lation would be expected to have a range of orbital periods. Both 
the uncertainty in the distances to UCWDs and in the formation 
kinematics restricts our ability to draw quantitative conclusions. 
However, both the observational samples presented in figures [2] 
&[3] with their different assumptions, seem to suggest that these 
apparently single LMWDs have experienced some kick. 

The broadly favoured period range at explosion is between 
10 and 100 days, with the longer orbital periods preferred for 
a larger initial scale-height. The sample presented by Vidrih et 
al. leads to a lower set of tangential velocities, and their sample 
strongly prefers 100 over 10 d. This may be because they assume 
fainter absolute magnitudes than we adopt, and hence systemati- 
cally produce lower velocities than us. Alternatively, their larger 
but less cool sample may be more likely to be contaminated by 
objects which are not single LMWDs. 

This period range of 10-100 d approximately encompasses 
the core masses appropriate for the production of LMWDs, ex- 
actly as might be expected if the stellar evolution was truncated 
by the loss of a giant star's envelope. This is consistent with 
an explanation of these UCWDs as being descended from giant 
donors. We note, however, that we cannot exclude a thick disc 
origin for these objects with an arbitrarily small kick. Although 
we cannot use this apparent kick as a definitive signature of an 
origin in type la supernovae, a thick disc origin would not falsify 
our SN la hypothesis, but it would suggest longer orbital periods 
at the point of explosion. 

5. Discussion 

The main conclusion of this study is that single LMWDs con- 
stitute indirect evidence that SNe la are formed through the 
single-degenerate channel, specifically from systems with red- 
giant donors. A field population of truly single LMWDs, espe- 
cially with masses less than ~ 0.4 M , would seem to require 
a SN la origin, or a significant revision of our understanding of 
either CV evolution or AIC03 That is based on evolutionary ar- 
guments (sections 2 & 3) and does not rely on any kinematic 
information. 

We have simultaneously demonstrated that there is a natural 
formation channel for single LMWDs which does not require 
any unexpected modification of single-star evolution. 

Those conclusions are independent of the nature of appar- 
ently single UCWDs. However the currently observed UCWD 
sample may be dominated by LMWDs, and we find that the kine- 
matics of the known single UCWDs seem to be consistent with 
an origin in SN la with red-giant donor stars. For reasonably 
broad assumptions, they even suggest that the orbital periods at 
the time of explosion were < 100 d. The population of UCWDs 
require further study in order to strengthen those specific con- 
clusions. 



Here 1— P x gives the probability that the data and model x are drawn 
from different distributions, so we can reject a model with P = 0.01 with 
99% confidence. 

17 It may be that explosions of sub-Chandrasekhar-mass WDs occur 
which are not seen as SN la but are capable of stripping the envelopes 
of RG donor stars. That would provide a further possible channel for 
the formation of single LMWDs. 
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5.1. SN la progenitor evolution 

The implication of the existence of single LMWDs is that at least 
some single-degenerate SNe la occur with the donor star on the 
giant branch at the time of the explosion. 

Both the high-velocity WD and single LMWD populations 
seem to contain members originating in SN la explosions. Dwarf 
star donors at the time of explosion should acquire higher space 
velocities than giant donors (due to their shorter orbital periods), 
but the less extended donors should be less stripped by the ram 
pressure of the explosion (Marietta et al., 2000) and hence go 
on to produce higher-mass WDs than a giant which loses its 
envelope prior to helium ignition. Hence the Oppenheimer et 
al. (2001) sample of high-velocity WDs considered by Hansen 
(2003) do not need to be low-mass for a SN la origin to be rea- 
sonable, nor do single LMWDs need to have high space veloci- 
ties to invite a SN la explanation for their production. LP 400-22 
(Kawka et al. 2006) is notable for being an extreme object using 
either selection criterion. 

A simplistic comparison of the remnant space densities 
quoted for the Oppenheimer et al. (2001) high-velocity sam- 
ple (1.8 x 10~ 4 pc~ 3 ) and the Gates et al. (2004) UCWDs (~ 
3 x 10~ 5 pc~ 3 ) suggests that the supersoft channel is almost an 
order of magnitude more important than the red-giant channel, 
though presumably some single LMWDs destined to be single 
UCWDs have not yet cooled sufficiently to become UCWDs, 
and some UCWDs may not be LMWDs. The relative impor- 
tance of these formation channels will be an interesting quan- 
tity to constrain with future data and compare with binary pop- 
ulation synthesis models. Perhaps more decisively, the estimates 
of Hansen (2003) and section 14.41 imply that the total number 
of observed remnants is consistent with the expected number of 
past SNe la in our Galaxy. If this is confirmed by future work, it 
would leave room for only a minority of SNe la to result from 
double-degenerate systems. 

One curiosity of our study of UCWDs is that the range of 
orbital periods we infer for the SN la progenitors at explosion is 
only broadly consistent with that predicted by Hachisu & Kato 
(2001) for systems with red-giant donors. Their models pro- 
duce no SNe la with final orbital periods 2.5 d < P or \, < 60 d; 
most of the final parameter space for red-giant donors has 
100 d < f or b < 1000 d. Our favoured assumptions seem to in- 
dicate orbital periods on the shorter side of 100 d, i.e. the kine- 
matic signal seems to be a little stronger than expected. A possi- 
ble resolution could be that our use of UCWDs is biased and 
preferentially selects systems with shorter orbital periods be- 
cause they produce lower-mass WDs, which cool more rapidly. 
Or perhaps the progenitor population is preferentially from the 
thick disc, and hence the remnants falsely appear kinematically 
hotter than in our models. It may also be that some non-LMWD 
single UCWDs from the thick disc or Galactic halo could be con- 
taminating these results. Alternatively, this may be an indication 
that the models need some modification; for example the sys- 
tems containing red-giant donors might produce SNe la earlier 
than expected, or perhaps a subset of systems from the supersoft 
channel does not explode until their donor stars have evolved 
more than the current models predict. 

5.2. The nature of UCWDs 

We encourage further work on the nature of UCWDs. It is im- 
portant to confirm that the apparently single UCWDs really are 
single. High signal-to-noise searches for any spectral features 
that allow radial-velocity measurements would be worthwhile. 



GAIA astrometry is a long-term hope for examining the single 
status of UCWDs, as well as providing accurate distances to all 
these objects. 

It is also important to understand the mass distribution of 
the UCWD population. The inferred mass of the best-studied 
UCWD, LHS 3250, (0.23 M [Bergeron & Leggett, 2002]) 
makes it a clear LMWD and hence a good candidate SN la 
remnant, but no other apparently single UCWD has a well- 
constrained mass estimate that we are aware of. 

5.3. Outlook: a definitive signature? 

A sample of non-UCWD single LMWDs has distinct advantage 
over the UCWDs in that the prescence of spectral lines allows 
for radial velocity measurements (see, e.g., Maxted et al. 2000). 
The significant difficulty for these objects is in selecting a large 
sample: whilst UCWDs stand out from survey photometry, hot- 
ter LMWD do not. Eisenstein et al. (2006) used the Sloan digital 
sky survey to produce a catalogue of over 9000 white dwarfs, 
and identify 13 WDs with masses < 0.3 M , of which 7 have 
masses < 0.2 M . These objects should be investigated for signs 
of a companion NS, or for radial velocity variations. 

With a substantial sample of single LMWDs, it would make 
sense to use directional information about the space velocities of 
the objects rather than just the magnitudes of the transverse ve- 
locities. We expect that the use of such information will help dis- 
tinguish between a large initial scale-height, long orbital period 
population and a thinner, shorter-period population. Our conclu- 
sions would be considerably stronger if we were sure about the 
initial kinematics of these stars; it is important to try and deter- 
mine whether the remnants come from the thick disc. 

Once the single LMWD sample becomes large enough, 
it could potentially be split into sub-samples with different 
WD masses. If our model is correct, then less massive single 
LMWDs should be kinematically hotter than more massive sin- 
gle LMWDs (see section l3~4l >. 

Since we began this work, van Leeuwen et al. (2007) and 
Kilic et al. (2007a) have searched for companions to LMWDs 
in the radio and optical wavebands, respectively. The search by 
van Leeuwen et al. of LMWDs for radio pulsations found none 
'down to flux densities of 0.6-0.8 mJy kpc 2 ', and concluded 
that 'a given low-mass helium-core white dwarf has a probability 
of < 0.18+0.05 of being in a binary with a radio pulsar'. For four 
WDs with masses < 0.2 M , Kilic et al. found: 'None of these 
white dwarfs show excess emission from a binary companion, 
and radial velocity searches will be necessary to constrain the 
nature of the unseen companions.' Our paper suggests that the 
assumption that there are unseen companions is not nessecary. 

5.4. Runaway hot subdwarfs 

Perhaps the most interesting possibility for the evolutionary state 
of the donor in a single-degenerate SN la is that of a hot sub- 
dwarf (sdO or sdB) star. This would naturally allow short or- 
bital periods (~ 1 h) and also naturally produce extremely low- 
mass WDs, as recently observed in the runaway WD LP 400-22 
(Kawka et al. 2006), which is inferred to have a mass of 0.17 
M and a tangential velocity of 414 + 43 km s . This would 
also provide a natural explanation for stars like the runaway hot 
subdwarf US 708 (discovered by Hirsch et al. 2005). We feel that 
a SN la origin for this object is more satisfying than a scenario 
combining dynamical ejection from the supermassive black hole 
in the Galactic centre with the simultaneous merger of two he- 
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Hum WDs (as speculated by Hirsch et al.). The orbital velocity 
in this case may well be augmented by a kick due to an impulse 
from the supernova shock (e.g. Marietta et al. 2000). 

The evolution of the WD-sdB binary system KPD 
1930+2752 (see, e.g., Maxted, Marsh & North 2000) has been 
investigated by Ergma, Fedorova & Yungelson (2001). They 
conclude that this system is likely to eventually result in a merger 
of two WDs (see also Geier et al. 2007). However we see no rea- 
son why similar systems could not produce a SN la via a single- 
degenerate channel, hence producing such objects as US 708 and 
then LP 400-22. 



6. Summary and Conclusions 

We have considered the formation of apparently single LMWDs 
in general, concluding that the most natural scenario for the for- 
mation of single LMWDs is that they are the remnants of donor 
stars in single-denenerate SNe la. Indeed, lone LMWDs should 
be expected if some single-degenerate SNe la do occur with gi- 
ant donor stars, as inferred from the observations of Patat et al. 
(2007), notably if the donors lose a significant fraction of their 
envelopes, as predicted for giant donors (Marietta et al. 2000). 
The observations of Maxted et al. (2000), van Leeuwen et al. 
(2006) and Kilic et al. (2006) are all in support of the existence 
of a population of genuinely single LMWDs. 

It seems difficult for the majority of apparently single 
UCWDs to posess companions, and we have adopted them as 
a useful sample of single LMWDs. We have integrated a popu- 
lation of SN la donor remnants through a simple Galactic poten- 
tial and compared the results of those calculations to the known 
space velocities of apparently single UCWDs. Our results are 
consistent with the single low-mass UCWDs having once been 
red-giant donor stars at the time of a SN la explosion, as pre- 
dicted for single LMWDs. 

A unified picture emerges in which the high-velocity WDs 
are remnants of main-sequence donors in SNe la (as suggested 
first by Hansen, 2003), and a kinematically cooler population 
of single LMWDs were once giant donors in long-period SN la 
progenitors: their longer orbital periods led to a lower runaway 
velocity wheras their tenuous envelopes were stripped more eas- 
ily by the supernova ejecta to produce LMWDs. 

Furthermore, it seems plausible that runaway LMWDs such 
as LP 400-22 and runaway hot subdwarf stars such as US 708 
originate from donor stars in short-period (~ 1 h) SN la systems. 
We will explore this idea in more detail in a future paper. 
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